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I. INTRODUCTION 

It is well known that the small- a; behavior of struc- 
ture functions of deep inelastic scattering is determined 
by the hard pomeron contribution. In the leading order 
the pomeron intercept is determined by the BFKL equa- 
tion [1] and the pomeron residue (the 7*7*-pomeron ver- 
tex) is given by the so-called impact factor. To find the 
small- a; structure functions in the next-to-leading order, 
one needs to know both the pomeron intercept and the 
impact factor. The NLO pomeron intercept was found 
many years ago [5] but the analytic expression for the 
NLO impact factor is obtained for the first time in the 
present paper. 

We calculate the NLO impact factor using the high- 
energy operator expansion of T-product of two vector 
currents in Wilson lines (see e.g the reviews p3H]). Let us 
recall the general logic of an operator expansion. In order 
to find a certain asymptotical behavior of an amplitude 
by OPE one should 

• Identify the relevant operators and factorize an am- 
plitude into a product of coefficient functions and 
matrix elements of these operators 

• Find the evolution equations of the operators with 
respect to the factorization scale 

• Solve these evolution equations 

• Convolute the solution with the initial conditions 
for the evolution and get the amplitude. 

Since we are interested in the small- X asymptotics of deep 
inelastic scattering (DIS) it is natural to factorize in ra- 
pidity: we introduce the rapidity divide r\ which separates 
the "fast" gluons from the "slow" ones. 

As a first step, we integrate over gluons with rapidities 
Y > r/ and leave the integration over Y < r\ for the later 
time, see Fig. 1. 
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FIG. 1. Rapidity factorization. The impact factors with 
Y > rj are given by diagrams in the shock-wave background. 
Wilson-line operators with Y < rj are denoted by dotted lines. 



It is convenient to use the background field formalism: 
we integrate over gluons with a > a = e n and leave glu- 
ons with a < a as a background field, to be integrated 
over later. Since the rapidities of the background glu- 
ons are very different from the rapidities of gluons in our 
Feynman diagrams, the background field can be taken in 
the form of a shock wave due to the Lorentz contraction. 
To derive the expression of a quark (or gluon) propagator 
in this shock-wave background we represent the propa- 
gator as a path integral over various trajectories, each of 
them weighed with the gauge factor Pexp(ig J dx^A^) or- 
dered along the propagation path. Now, since the shock 
wave is very thin, quarks (or gluons) do not have time to 
deviate in transverse direction so their trajectory inside 
the shock wave can be approximated by a segment of the 
straight line. Moreover, since there is no external field 
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outside the shock wave, the integral over the segment 
of straight line can be formally extended to ±00 limits 
yielding the Wilson-line gauge factor 
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du p^A a (up\ + x±) 
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where the Sudakov variable otk is defined as usual, k = 
ctkPi + PkP2 + k± . We define the light-like vectors pi and 

P2 such that q = pi — x bPi and p — P2 + ^-pi where q 
is the virtual photon momentum, p is the momentum of 
the target particle, and xb — Q 2 /s <C 1 is the Bjorken 
variable (at large energies s ~ 2p ■ q). The structure of 
the propagator in a shock-wave background looks as fol- 
lows (see Fig. [2]): 

[Free propagation from initial point x to the point of in- 
tersection with the shock wave z\ 

x [interaction with the shock wave described by the 
Wilson-line operator C/ z ] 

x [Free propagation from point of interaction z to the 
final point y\ . 



I U(z ± ) 



FIG. 2. Propagator in a shock-wave background 




FIG. 3. Impact factor in the leading order. Solid lines repre- 
sent quarks. 
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Here we introduced the conformal vectors [5J |S] 
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and the notation 1Z 



« 2 (Cl'C2) 



The above equation 



2(«-Cl)(K-C2)' 

is explicitly Mobius invariant. In addition, it is easy to 
check that ^^-(r.h.sWO. 

Our goal is the NLO contribution to the r.h.s. of Eq. 
Q, but first let us briefly discuss the three remaining 
steps of the high-energy OPE. The evolution equation 
for color dipoles has the form [SJ [TO] 



The explicit form of quark propagator in a shock- wave 
background can be taken from Ref. [5] 
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As usual, we label operators by hats and (O) a means 
the vacuum average of the operator O in the presence 
of an external field A. Hereafter use the notations — 



fX 

P2 x fJ- 



_ Vi, 



- PlXfi 



(and our metric is 

(1,-1,-1,-1)). Note that the Regge limit in the coordinate 
space can be achieved by rescaling 

2 2 

X -> px* -pi + X, — £>2 + £_L, 

s sp 
2 2 

v -> py*-pi + y.—P2 + y±, (3) 

s sp 

with p — y 00 , see the discussion in Refs. [BJEI- 

The result of the integration over gluons with rapidities 
Y > 77 gives the impact factor - the amplitude of the 
transition of virtual photon in two- Wilson-lines operators 
(sometimes called "color dipole" ) . The LO impact factor 
is a product of two propagators |2]), see Fig. [3] 
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(To save space, hereafter Zi stand for Zi± so zf 2 = z 
etc.) The explicit form of the NLO contributions can 
be found in Refs. [U [H] Q2] while the agrument of the 
coupling constant in the above equation (following from 
the NLO calculations) is discussed in Refs. ( |13l 114] ). 

It is worth noting that we performed the OPE pro- 
gram outlined above for scattering of scalar "particles" 
in J\f = 4 SYM and obtained the explicit expression for 
the four-point correlator of scalar operators at high ener- 
gies in the next-to-leading order [7] . In QCD the analytic 
solution of the evolution equation for color dipoles with 
running coupling constant is not known at present. This 
prevents us from getting the explicit NLO amplitude as in 
Af = 4 case. We can, however, perform the first two steps 
in our OPE program discussed in the Introduction: calcu- 
late the coefficient function (impact factor) and find the 
evolution equation for color dipoles. The next two steps, 
solution of the evolution equation ^ with appropriate 
initial conditions and the eventual comparison with ex- 
perimental DIS data are discussed in many papers (see 
e.g. [TS]). It is worth noting that, contrary to the evolu- 
tion equation, the NLO correction to the impact factor 
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has nothing to do with running of the coupling constant 
- it starts at the NNLO level. Thus, the argument of the 
coupling constant at the NLO level is determined solely 
by the evolution equation for color dipoles. For numer- 
ical estimates involving the impact factor one can take 
a s (\x — y\) as the first approximation since the char- 
acteristic transverse distances in the impact factor are 
~ \x- y\. 

The paper is organized as follows: in Sect. 2 and 3 we 
calculate the NLO impact factor in the coordinate rep- 
resentation (the results of these Sections were published 
previously in Brief Report [E]). The Mellin representa- 
tion of the impact factor is presented in Sect. 4 and Sect. 
5 contains the impact factor in the momentum represen- 
tation for the forward case corresponding to deep inelas- 
tic scattering. Finally, we present the NLO BFKL kernel 
and discuss the fc^-factorization for DIS in Sect. 6. 





FIG. 4. Impact factor in the next-to-leading order. 



currents in Af = 4 SYM carried out in our previous pa- 
per [12] . The gluon propagator in the shock- wave back- 



ground at > > in the light-like gauge 
is given bv [171 ITS] 







II. CALCULATION OF THE NLO IMPACT 
FACTOR 



Now we would like to repeat the steps of operator ex- 
pansion discussed above to the NLO accuracy. A general 
form of the expansion of T-product of the electromag- 
netic currents in color dipoles looks as follows: 

(x - yfT{Ux)l^{x)Uy)l v ^{v)} 
' d 2 z\d 2 z 2 
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{l™( Zl ,z 2 )[l + ^]tr{UlU%} 



x MUlUl:}tv{UlU^} - NMUlUl:}]} (7) 

Unfortunately, in terms of Wilson-line approach there is 
no direct way to get the NLO impact factor for the BFKL 
pomeron. One needs first to find the coefficient in front 
of the four- Wilson-line operator (which we will also call 
the NLO impact factor) and then linearize it. 

The structure of the NLO contribution is clear from 
the topology of diagrams in the shock-wave background, 
see Fig. |4]below. Also, the term ~ 1 + can be restored 
from the requirement that at U = 1 (no shock wave) one 
should get the perturbative series for the polarization 
operator 1 + ^ + 0(a 2 s ). 

In our notations 



l)iv z i) 
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t 6 (k-Ci)(k-C 2 ) 
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x ^i;[(-Ci)(-C 2 )™2 k2 ^-^]- ( 8 ) 

which corresponds to the well-known expression for the 
LO impact factor in the momentum space. 

The NLO impact factor is given by the diagrams shown 
in Fig. [4] The calculation of these diagrams is similar 
to the calculation of the NLO impact factor for scalar 
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1^ ~ P2»(X - Z)f ^ 1 yj^-p 2v (y-z) 



ir 2 [{x — z) 2 + it 



z± q(z) ir 2 [(z - y) 2 + it] 2 



where -k- can be either a -. . or a 1 . 



which leads to 



the same result. (This is obvious for the leading order 
and correct in NLO after subtraction of the leading-order 
contribution, see Eq. (16) below). 



The diagrams in Fig. [4] can be calculated using the 
conformal integral 



d z — 7^- 



(x-z) 4 ^{z-y) 4 z 4 



H o v 



x 2 y 2 (x — y) 2 



- fi <-> v (10) 



which gives the 3-point ipipF^ Green function in the 
leading order in g. Using Eqs. ([9| and (10), per- 

forming integrals over z.'s and taking traces one gets 
after some algebra the NLO contribution of diagrams in 
Fig. [4] in the form 

i F ^ Zl ,z 2 , Z3 ) = rr(zuZ2,z 3 ) + r^( Zl ,z 2 ,z 3 ) (11) 

where 



I\ {zi,z 2 ,z 3 ) 
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and 

{h)iiv{Z\,Z 2 ,Z 3 ) = 

x / (« • C2) d 2 
1 (k • C3) dx^dy" 
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(recall that z 2 




2(C 2 • Ca) 



'ij± ~ 2(Ci'Ci) andZ, = ^(k-CO). We 
obtained this expression at 1, > > but from the 
conformal structure of the result it is clear that this ex- 
pression holds true at a;* < < as well. 



The integral over a in the r.h.s. of Eq. ( 12 ) diverges. 
This divergence reflects the fact that the contributions 
of the diagrams in Fig. [4] is not exactly the NLO impact 
factor since we must subtract the matrix element of the 
leading-order contribution. Indeed, the NLO impact fac- 
tor is a coefficient function defined according to Eq. Q. 
To find the NLO impact factor, we consider the opera- 
tor equation ^ in the shock-wave background (in the 
leading order (U Z3 )a — U Z3 ): 

(T{^(xh^(x)kyh u ip(y)})A 

d 2 Z\d 2 Z 2 j2 r NLO I \ 

a « z 3 V {x,y;zi,z 2 ,z 3 ;r]) 



(14) 



[tr{U Zl Ul }tr{U Z3 Ul } - N c tr{U Zl 17*, }] 



The NLO matrix element {T{'4){x)^^{x)ip{y)Y^{y)}) A 
is given by Eq. ( 11 ) while the subtracted term is 
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'0 " J ^13^23 

x MU Zl Ut 3 }tr{U Z3 Ul} - NMUzMj] (15) 

as follows from Eq. ([6]). The a integration is cut from 
above by a — e ri in accordance with th e de finition of 
operators U n , see Eq. (|l|. Subtracting (15 1 from Eq. 



( 11 1 we get 
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where C is the Euler constant. Note that one should ex- 
pect the NLO impact factor to be conformally invariant 
since it is determined by tree diagrams in Fig. H| How- 
ever, as discussed in Refs. [HEHH], formally the light- 
like Wilson lines are conformally (Mobius) invariant but 



the longitudinal cutoff a < a in Eq. (1) violates this 
property so the term ~ In aZ 3 in the r.h.s. of Eq. ( 16 1 is 



not invariant. As was demonstrated in these papers, one 
can define a composite operator in the form 



Mu Zl ut 2 }] a 
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where a is an arbitrary constant. It is convenient 
to choose the rapidity-dependent constant a 



ae 



-2*1 



so that the [tr{f/^ 1 C/^}]™ nf does not depend on 
r) = ln<7 and all the rapidity dependence is encoded 
into a-dependence. Indeed, it is easy to see that 

iMU Zl Ut 2 }] C a ° nf = and A M U Zl Ul 2 }] c : n! * dc " 
termined by the NLO BK kernel which is a sum of the 
conformal part and the running-coupling part with our 
0(a 2 ) accuracy [HH]. 



Rewritten in terms of composite dipoles (17 1, the op- 
erator expansion Q takes the form: 

T{^(x)rHx)kyh^(y)} 

'd 2 z 1 d 2 z 2 r LO 
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We need to choose the "new rapidity cutoff" a in such 
a way that all the energy dependence is included in the 
matrix element (s) of Wilson-line operators so the impact 
factor should not depend on energy. A suitable choice of 
a is given by ao — —n~ 2 + ie 



s(x-y) 2 

(x - v fT{Ux)i^{x)Uy)i^{y)} 

' d 2 z\d 2 z 2 
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Here the composite dipole [tr{?7f U]£}]a a is given by 
Eq. Q with a = -j^jj + it while I^( Zl ,z 2 ) and 
I 2 "(zi,z 2l z 3 ) are given by Eqs. ^ and (131, respec- 
tively. 



III. NLO IMPACT FACTOR FOR THE BFKL 
POMERON 

For the studies of DIS with the linear NLO BFKL 
equation (up to two-gluon accuracy) we need the lin- 
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earized version of Eq. (19). If we define 
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and consider the linearization 
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one of the integrals over Zi in the r.h.s. of Eq. ( 19 ) can 
be performed. The result is 
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where Li 2 (z) is the dilogarithm. Here one easily recog- 
nizes five conformal tensor structures discussed in Ref. 

While it is easy to see that 



1 



llf 3 {x,y:,zi) = (24) 



dx„ (x - yf dx» dy v aP 

one should be careful when checking the electromagnetic 
gauge invariance in the next-to-leading order. The reason 
is that the composite dipole U a °(zi, z 2 ) depends on x via 



the rapidity cutoff ao = — ,f x ^ v \ 2 so from Eq. (21) we 
get 



^Jx^Wdx^WJ ° { ' } aP {,V,) 



= 2 
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Using the leading-order BFKL equation in the dipole 
form (linearization of Eq. ([6])) 
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x [U a (z 1 ,z 3 ) +U a (z 2 ,z 3 ) -U a {z 1 ,z 2 )] (26) 
we obtain the following consequence of gauge invariance 
d 1 dn a dn 
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We have verified that the expression (23) satisfies the 
above equation. 



IV. PHOTON IMPACT FACTOR IN THE 
MELLIN REPRESENTATION 



In preparation for Fourier transformation we calcu- 
lated the Mellin transform of the photon impact factor 
(23). We project the impact factor on the conformal 



eigenf unctions of the BFKL equation [20] 
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(here z — z x + iz y ,z — z x — iz y , z w = Z\ — z etc.). 
Since electromagnetic currents are vectors, the only 
non- vanishing contribution comes from projection on 
the eigenfunctions with spin and spin 2. The spin-0 
projection has the form (throughout the paper we 
reserve the notation 7 for J + iv): 
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where 7 = i + ii/, 7 = 1 — 7=5 
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Using the decomposition of the product of the trans- 
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The contribution of spin 2 in the t-channel has the form is a composite dipolc ( 17 1 in the Mellin representation. 
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get the high-energy OPE in the form 



Substituting the decomposition (35)) in Eq. (21) we 

(38) 



1 



{x-yfT{f{x)f{y)} 



% J d2z o p [(! + y\ y) 

+ (^ 2 + l)[{(l + ^)j 2 %(*,y;zo^) 

+ J2,a f (X, V, Zq, v) }U aa (V, Zq) + 



The Eq. ( 38 ) and its Fourier transform ( 46 ) are the main 



results of this paper. 

At this point it is instructive to check again the photon 
gauge invariance ^T{j^{x)j u {y)} = 0. Since a = 
— k 2 + ie we need to differentiate U aa (y, zq) too: 



d - d - 

—U ao (v,z ) = a—U a (p,z ) 



d 
dx^ 



lna (39) 



Let us start with spin-2 contribution. It is easy to see 
that 



d S£ v 
dx^ (x — y) 6 



0, 



d 

S ^dx^ lna ° = ° 



and therefore the second term in r.h.s. of Eq. (38) is 
gauge invariant (recall that ao = — kT 2 + ie). 

For spin-0 part we need to use eq. ( 39 1 . Since U a (y, zq) 
are the projections of color dipoles on the eigenfunctions 
( 28 1 of the BFKL equation, the evolution equation ( 26 ) 



simplifies to 



2a—U a (v 1 z ) = uj{v)U a {v,z Q ) 
da 



(40) 



where u){v) — as ^ c x 7 is the BFKL pomeron intercept 
(as usual 7=5+ w). We obtain 

—U ag (v,z ) = -^U ao (v,z )—\n K 2 (41) 



In the leading order the derivative ( 41 ) does not con- 



tribute so the formula for gauge invariance is simply 







J^(x,y;z ,u) = 



dx^ 

It is easy to demonstrate that J^Q in the r.h.s. of Eq. 



(29) satisfies this requirement. 



In the NLO we need both and u J^9 parts so the 

requirement for electromagnetic gauge invariance takes 
the form 

d _h[ n / \ , n . , <91iik 2 



(42) 

We have checked that the r.h.s. of Eq. ( |29[ ) satisfies this 
equation. 



PHOTON IMPACT FACTOR IN THE 
MOMENTUM SPACE 



In general, the rapidity evolution of color dipoles is 
non-linear but in this paper we assume that we can lin- 
earize it to the dipole form of the BFKL equation, like 
in the case of scattering of two virtual photons. More- 
over, we will consider only the forward case which cor- 
responds to deep inelastic scattering. In this case, one 
may write down the high-energy OPE in the form of fcy- 
factorization formula 



where 



d 4 x e^{p\T{j^x)U0)}\p) 



d 2 fc±V(g,fc ± )((p|W(fcL)|p)) (43) 



(U(k ± )) = dx ± e- i ^ x ^(U(x ± ,0)), 



q = Pi + '^rP2 and p = P2 + —Pi is the target's momen- 
tum. The reduced matrix element ((p\U(k)\p)) is defined 



(p\U(h)\p + Pp 2 ) = 2 V S(P)({p\U(k)\p}) 
((p\U(k)\p)) = fd 2 ze^ k ^((p\U(z,0)\p)) (44) 



where the factor 2iv6(/3) reflects the fact that the for- 
ward matrix element of the operator U X U^ contains an 
unrestricted integration along the p\. Our goal in this 
Section is to find the impact factor Ip, v (q, k±) in the next- 
to-leading order. 

Since our "energy scale" ao = —k~ 2 for color dipoles 
depends on x an d v . to perform the Fourier transforma- 



tion of the OPE ( 38 1 one should express U ao in terms of 
Ua m with a m independent of coordinates x and y. A suit- 
able choice is a m — 1/xb- With this choice, the impact 
factor does not scale with s and all the energy dependence 
is included in martix elements of color dipoles. This is 
similar to the choice /i 2 = Q 2 for the DGLAP evolution: 
the coefficient functions in front of the light-ray opera- 
tors will not depend on Q 2 (except for a s (Q 2 ) of course) 
and all the Q 2 dependence is shifted to parton densities. 
The leading-order evolution of a color dipole U a is given 
by Eq. p| 



Ua {v,ZQ) = U am (v,z Q )(aox B ) 2 , 
U^{v,z ) = ^ 2 >,z )(a x B )^\ (45) 



so the Fourier transform of Eq. ( 38 ) yields 



1 



d xd y 5(y,) e 



iq-(x-y) 



T{^^{x)ip^ y ip{y)} (46) 



dv 

27-1 



d Zq 



r r (7 



u(v) \ r 2 



r 2 (2- 7 



lr(4-27 + w(i/))r(2 + 7 +^) 



20 

7r(3-7)r( 7 



r(2- 7 ) ( 7 7 + 2)Pr 1 + - + 

1 L V 7T 



2tt 



+ (3 7 7 + 2)Pr(l + ^ + ~^{v))]u an {z 0} v) 

^)r 2 (2- 7 + ^) 



2r(4 - 2 7 + w(2, i/))r(2 + 7 + ^) 

ri(2) _ i(2) 

W 0n >,,*o)P'" / + U a (v,z )P" 1 



X 1 + 



a, TV, 



2tt 



J r(2 7 ) 



■■^ 6( „ )n ^mi±i) 



we get 



d 4 ^* e iqx (p\T{j fl (x+-z*p 1 )j v (-z* Pl )}\ P ) 



N, 



dv 



^ (off) 



fc 2 ^r 3 (7)r 3 ( 7 ) (p|w(fe)| P ) 



{ 



27 32tt 3 7 v q 2; r(2 7 )r(2 7 ) 3 + 477 
r( 7 + ^M)r 2 (2 - 7 + ^M)r(4 - 2 7 )r(2 + 7) 
r( 7 )r 2 (2 - 7 )r(4 - 2 7 + w(i/))r(2 + 7 + ^) 

(77 + 2)Pr(l + ^ + ^*iM) 

+ (377 + 2)Prfl + — + ^$2M 

V IT 



2tt 



r( 7 + ^)r 2 (2- 7 - 



w(2,i/) 
2 



)r(4-2 7 )r(2 + 7 ) 



r( 7 )r 2 (2 - 7 )r(4 - 2 7 + «(2, i/))r(2 + 7 + 



77 



7T 



a s JV c 
2tt 



(1/))} (50) 



where 



where 







q 2 


p r 




A 1 2 1 
P2<T \ 








= Of 1 - 






= (.^ + 





p 2 v 



9 -P2 



ig" 



and 



F( 1 ) + ^ ^_ 

w; 2 + 77 18(2- 7 ) 



1 



1 



7 



10 



27 2 7 18(1+7) 3(1 + 7 ) 2 

F(l) + ±- 
277 

7 X 7 , X 7 ( X + 3 7) 

2(2 + 377) 1 + 7 2 + 377 



(47) 



(48) 



The last step is to rewrite Eq. ( 46 ) in the /cT-factorized 



form (43). Using 



[k 2 P^ + fc 2 P^] - 2 K± x 



fei 



There is a subtle point in in the Fourier transform of Eq. 
( 46 1 : the contribution of infinite zq doubles the result 



from finite zq as shown in the Appendix. 

Note that w(u) = ^x 7 and w(2, v) = ^(x 7 



77- 



are of order a s so one should expand the expressions 



like T(4 — 27 + uj{v)) in Eq. (46) up to the first order 



in lo{v). Using definitions (43) and (44) one obtains the 
impact factor in the form 



N r [dv 



sinh ', 



1287 ^ (1 + ^ 2 ) cosh 2 



("+- 2 ) 



7T 

o „ 



a,N c 
2tt 

a s N c 
2tt 



\q 2 ) 

Fx(v)]pr 



r 2 



1 



V 

~2 



7T 



27T 



(u)]pf) (51) 



where (as usual, 7 = | + «V) 

■^l(2)M = *l(2)(» / )+X7*(« / )i 

,(1/) = F 6 (i/) + (x 7 -^-)*M 



(52) 



^3 



= ^(7) +2^(2 -7) -2^(4 -2 7 ) -^(2 + 7) 



v 1 °' 7/ 4^r(2 -7)r 2 (7)r(2 7 )7 v v 7/ ' 
( V*)> = - r(7)r ( 3-27)r 2 (i + 7) 



47r(3-7)r 2 (2- 7 )r(i + 2 7 ) 

£ - (2) 

x I d 2 k ^{u (k)) 



(49) 



The structures P\ and P2 correspond to unpolarized 
structure functions F\(xb) and i^C^s)- The third term 
vanishes for nucleon structure function but contributes 
to polarized structure functions of a vector meson (or 
photon) . 



It is instructive to compare Eq. (51) with the well- 



known double-integral representation of the leading- 



9 



order impact factor (see e.g. Ref. [ST], [3]): 



I^(q,k) 



32 J Q 2 uu + kj_vv 



dudv- 



{< 



1 — 2uu 



— 2vv + 4muvv)P^ v + (1 — 2uu — 2vv + 12uuvv) P£ 

+ 2™^- 2 + ^~^| (53) 

(Note that the definition of I^(q,k) differs in sign and 
N c from that of Ref. [5] ) . It is easy to see that Eq. ( 53 ) 
is equal to LO terms in the r.h.s. of Eq. (51 1. 



With the fcr-factorization in view let us rewrite the 



evolution equation ( 54 ) in terms of 



Va m (z±) = -d 2 x {{U am {zx,0)}) (57) 

proportional to the dipole unintegrated gluon distribu- 
tion T>{xb, z_l, fi) 

47r 2 X£> 

V XB (z±,fi) = — — a s (fJ,)V(x B ,z ± ,n) (58) 

where 



VI. NLO BFKL FOR COLOR DIPOLES 

For completeness, in this Section we present (lin- 
earized) evolution equation for composite color dipoles 
and discuss how it is related to usual NLO BFKL ap- 
proach [5j . The evolution equation for forward matrix el- 
ements of color dipoles U(z 12 ) = (W(zi, z 2 )) reads [T11IT2"] 



„ d , , , . a s N c 



d 2 z' i 1 



— f] 

,/ 2 \ 



+ 



4?r 



z'\z~ z >y 

V A2 „/2 



,67 tt 2 , 1Q 



}[u a (z')+U a (z-z')-U a (z) 



+ 



47T 3 



j2 ' 

a z 



1 



In 2 



(z-z') 2 z' 



+ Fi -. :' i + <!>( :. :') W a (/) + 3^((3)W a (z) (54) 



where 



F(z,z') = 1 



NV I6z 2 z' 2 



" 2 -2z 2 z' 2 / 2 



z 2 z' 2 



z 2 -z' 2 ^ z 2 
m — k 



z 2 z' 2 



1 



3z 4 + 3z' 4 -2z 2 z' 2 



'\2 



16z 4 z 



4j4 



(z,z') 



N 3 

oo 



1 - 

dt 



{z 2 + z' z ) 



iz 2 z /2 



In- 



1 + t 



and 

<D(z,z') = 



(z 2 - z' 2 ) 



o z 2 + i 2 z' 2 

(55) 

2 z 2 z' 2 (z-z') 4 



In — ?r In 



(z-z') 2 (z + z') 2 
+ 2Li 2 (-^)-2Li 2 (-^)]- ( 



(z 2 + z' 2 ) 



2n 4 



1 



(z 2 



(z-z') 2 (z + z') 2 



Z' 

1 /"OO 



(56) 



(z — z'w) 2 



In 



2 /2 

U Z 



Note that the kernel is a sum of the "running-coupling" 
part proportional to b — ^N c — |nj and the conformal 
part, see the discussion in Ref. [1]. Here a s = a s (/i) and 
(i is the normalization point in the MS scheme. 



T>(x B , zx,fi) 

= ( 2 ) 2 /^^(plTrj [oopi + zj_, 2 z*pi + zj_l 

S J TtXB S J 

2 2 
F.f (~2*|>i + zj_) [-a;*pi + Zj_,-oopi + zj 

x [-oppi.O^.^Oj^oppi]}!!?)"™^ 

Hereafter we use the notation 



(59) 



[x,y] = Pexp^igJ du (x - yY A^ux + uy)j 

for the gauge link connecting points x and y. The color 
dipole is renorm-invariant so T> depends on \x to compen- 
sate g 2 {p) dependence. 
The Fourier transform 

V{x B ,k±) = J dzxe l{k ^V(x Bl zx) 

is called the dipole gluon TMD (transverse momentum 
dependent distribution). Note, however, that the dipole 
gluon TMD defined above differs from the definition 



T>(x B ,zx,f-i,v) 
,2,2 f dz* 

J irx B 

2 2 
F.e(~z*px + z ± ) + zj_,-oopi + z ± ] 



-oop 1 ,0]F.e(p)[0,oop 1 ]}\p)i 



(60) 



which reduces to the usual parton density at zx = 0. It 
should be emphasized that Eq. ( 60 ) is a more complex 



operator than (59). The difference is especially clear in 



the case of Af = 4 theory: the dipole gluon TMD ( 59 ) 



is UV finite while Eq. (60) is UV divergent so it needs 



additional UV counterterms, see the discussion in |22) . 
(These UV divergent terms are directly proportional to 
xb so they vanish for the definition (59)). Also, the role 
of parameter xb is different in the two definitions: in Eq. 
([59l it is defined as a rapidity cutoff a m while in Eq. ([60l 



the rapidity cutoff r\ should be imposed separately from 
Xb- 



Differentiating Eq. ( 54 ) two times with respect to z 



we obtain the NLO BFKL evolution for dipole gluon 
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TMD ( 59 ) in the form 



2a — V a (z) 
da 

_ a s N c f 2 



«'V<{(l 



a s b 

47T 



4 96 



n 2 N c 10n/-|\ r V a (z') (zz')V a (z) 
36 96 J/ L (z-z') 2 z' 2 (z-z') 2 

Os6 Va(zQ-Va(z) (z-z') 2 

2tt (z-z') 2 11 ' 



v/2 



a s 7V c r ln 2 (z 2 /z' 2 ) 
(z-z>) 



4tt 
a 2 7V 2 

3 i^rC(3)v a (z) 



^+ F(z,z') + $(z,z')]v a (z')} 

(61) 



Next we need to perform the Fourier transformation of 



Eq. (61). It can be demonstrated that 

InV/g' 2 ) 



d 2 qd 2 q' <(g )Z )-i( 9 ', z ') 
4tt 2 



+ F{q,q') 



+ <%></)] = -^^P+F(z,z') + $(z,z')(62) 

[z — z 1 y 

so the conformal part of the kernel looks the same in 
coordinate and momentum representations. 

Performing also the Fourier transformation of the 
running-coupling part one obtains the momentum- 
representation kernel in the form 



2a— V a (fc) 

da tt 



Vg(k') (k,k')V a {k) 

(k-k') 2 fc' 2 (fc-fc') 2 



X 1 + 



4tt 



In 



/i 2 N C( &7 it 2 10n/, 

k 2 + ~6~'"9~ ~ T ~ ~m~ r 



>]) 



Air 



V a (k') (k-k'f k 2 V a (k) (k-k') 2 

m k s In 



(k-k')" 2 



k' z (k-k') 2 



k 2 



4tt 
, a 2 s N 2 
2tt 2 



ln 2 (fc 2 /fc' 2 ) 
(k - k>) 



*U + F(k,k') + ^(k,k')] V a {k')) 



C(3)Va(*) 



(63) 



where V(Jfe) = Jdz e~ i{k ' z ^V(z). 



In terms of Mellin projections (7 = h + iv as usual) 



V W = £/ 2^V(n,i/)(* 2 r- 1 (A/A) 

n=CT 

V(n,i/) = fd 2 k (k 2 )-*(k/k) n/2 V{k) 



n/2 



r(7 + f) 



the kernel ( 63 1 takes the form 



2a— V a (n, i/) 
aa 



(65) 



a s l\ c r a s iV c r b / 2 

+ /(n, 7 )]}v o (n,j/) 



f/7 



X("»7) , x'fo,7) 
2 2 X (n >7 ) 



Here 

/(»»7) = [y ~ y ~ yyM n >^) - *"( n >7) 
+ F(n, 7) - 2$(n, 7) - 2$(n, 1 - 7) + 6C(3) (66) 

where 



'On 



77 



COS 7T7 



iY c 3 ^2(3-27)(l + 2 7 )""' 1 J (1 -27)sin 2 7r 7 



and 



$(n,7) 



7 1 + t 1 12 2 V V 



7r 2 1 , , / n + 1 



-Li 



1 + * 

(t) - Li 2 (-t) - [^(n + 1) - V(l) + ln(l + t) 
(-t) k 



*;=i fc=i v ; 

To compare to NLO BFKL from Ref. [2] one should 
rewrite above equation in terms of 



C(k) = ^vV(fc) 
g 2 (k) 



(68) 



since two gluons in the dipole U(k) come with extra g 2 
factor. The Eq. ( 65 ) turns into 



2a — C a (k) = 
da 

(k,k')C a (k) 



a s (k 2 )N c 



d 2 k 



Cg{k>) 

[k-k') 2 



k' z (k-k') 2 

C a (k') 



4-7T 

a s N c 

An 
„a 2 N 2 



(k-k') 2 k' 2 (k-k') 2 



a s N c /67 tt 2 10n/\ 
4tt V~9~ ~ ~3~ ~ 9N C ) 

k 2 C a (k) ! (k-k'f 

m 



k 2 



ln 2 (fc 2 /fc' 2 ) 
(fc - fc') 



— 2+ F(k,k') + <S>(k,k')] C a (k')) 



C(3)A(fc) 



2tt 2 

with the eigenvalues 



(69) 
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a s {k 2 )N c . , 



d 2 k'[ 



(k' 2 /k 2 )^e m t 



(k,k>) 



{k-k') 2 k' 2 (k-k') 2 



a s N c ,67 ir 2 10n fs 

1 H ! — H -) 

4tt v 9 3 9N C ' 



+ 



An 

a s N, 
An 



(k' 2 /k 2 )-^e ln ' 1 ' 



{k-k') 2 k' 2 (k-k') 2 
ln 2 (fc 2 /fc' 2 ) 



In 



(k - k') 



'\2 



k 2 



(k - k') 2 



F(k,k') + ®{k,k'] 



(fc' 2 /k 2 )~*>e m *X + 3%^C(3) 



a s (k 2 )N, 



°-{x{n,l) 



+ 



a s N c 
An 



b 

2N, 



[ X 2 (™,7) + x'(™,l)]+/K7)]} (70) 



which coincide with eigenvalues of the kernel [5] of the 
partial wave of the forward reggeized gluon scattering 
amplitude 

u J GM,q / ) = s {2) (q-q')+ [d 2 P K(q 7 p)GUp,q r ) (n) 



2 ^ 

a.,N, 



a s (q) 



+ < ^{lM-—[x\n n )+x 2 {n^)])] (72) 

This is somewhat surprising since the evolution of the 
composite (in Af = A SYM - conformal) dipole with re- 
spect to a gives the evolution of forward reggeized gluon 
scattering amplitude with respect to rapidity rj (of which 
<jj is the Mellin transform). To illustrate the transition 
between the two evolutions let us consider the calculation 
of the dipole evolution directly from the NLO BFKL for 
reggeized gluons. 

The impact factor &A(q) for the color dipole U(x, y) is 
proportional to a s (q)(e tqx — e lqv )(e~ lqx — e~ lqy ) so one 
obtains the cross section of the scattering of color dipole 
in the form 



W{x) = 



d-q d q' a s (q) 



q 2 q' 2 An 2 



(e iqx - l){e~ iqx - l)$ B {q') 
x / Pr(^r) G u (q,q') (73) 

QbW) is the target impact factor. To get the evolution 
equation with respect to rapidity one should change the 
energy scale to q' 2 . 



W(x) 



d 2 q d 2 q' a s (q) 



q 2 q' 2 An 2 



/>a+ioo 



d'jj 



2ni V q' z J 



(e lqx - l){e- iqx - 

(74) 



where G u (q 1 q') is the modified kernel with eigenvalues 

shifted by 2 X {n,v) X ' {n,v)^£^ , see Ref. [2]. The corre- 
sponding equation for V(n, v) takes the form [TT] 



d a s N, 
— V'{n,v) = 
drj 



n I 



ln^ 2 



d 
d-< 



a s N t 
An 

X(n,l) , x'( n ,l) 



N, 



2 2x(n, 7 ) 
+ f(n, 7) + 2 X (n, u)x\n 7 v)] }v"(n, v) 



(75) 



Let us demonstrate that it agrees with Eq. (65). The 
Mellin projection of composite dipole can be obtained 
from Eq. (Tl7|): 



V a (n,v) = (1 + 



a K N n 



An 



X 2 {n,i) + 3x'(n,7) 



+2x(n, 7 )(lna-2r; + 2C- 



dj' 



r a,N c 



v An 
+ 2(lna - 277) 
X'(n,j) 



) 2 {2 X 2 (n, 7 )(lna-27? 
b 



-X(n, 7 )(ln^ + ^ 



X{n,i) 



+ /(n,7) + 3x(n,7)x'(n,7) 



2 x(n,j)< 
+X 3 (n,7)] +X(n,v)}w(n,v) 



(76) 



The first-order term can be derived from Eq. (17) 
while the term ~ a 2 can be restored from the condition 
^V a (n, v) = up to an unknown function X(n, v) which 
requires NNLO calculation (and does not contribute to 
the NLO evolution). Now one can see that the deriva- 
tive with respect to a gives Eq. (65). Thus, the transi- 



tion between evolution of the composite dipole V a {n, v) 
with respect to a and the rapidity evolution of the dipole 
V (n, v) corresponds to the shift in eigenvalues on the 

function 2x(n, v)x'(n, v) 



the same transition that 



describes the shift of eigen valu es when going from energy 



scale qq' to q' in formula ( 74 ) 



VII. CONCLUSIONS 

Let us present again the kx factorization formula for 
DIS in the next-to-leading order: 



d*x e* qx {p\T{ Jtl (x)UmP) 



d 2 kj_ 



Inv{q, k±)Va m =x B (kj_) 



(77) 



where I^(q,k±) is given by Eq. (51) and the evolution 
equation for V a (k±) by Eq. (63) (or Eq. (651 in the 
Mellin representation) . 

The analytic NLO photon impact factor in momentum 



space for the pomeron contribution (51 ) and the NLO fcx 



12 



factorization formula ( 
ing are the main resu 



77) for the deep inelastic scatter- 
Previously, the 



ts of this paper 
impact factor was known only as a combination of an- 
alytic and numerical results [53]. It should be noted, 



Consider now the Fourier transform of this equation for 
the case of forward scattering 



however, that our final NLO result (51) is defined as a 



coefficient function in front of a composite operator (17) 



rather than in front of a usual dipole (with the rapidity 
cutoff) which corresponds to the impact factor defined in 
Ref. [2S] papers. 



f d A xd A y 8{y.) e iq < x - y) 
(x-y) li (x-y) 



9^ 



(x - yf 



dz\dz 2 TZ 2 

42 ( x - v) 6 

U{z 12 ) 



(80) 



It would be also instructive to compare our result ( 19 ) 



for the coefficient in front of the four- Wilson-line opera- 
tor (relevant for the structure functions of DIS off a large 
nucleus) to similar result for the NLO impact factor ob- 
tained recently in Ref. [24] using the dipole model. How- 
ever, as we already mentioned our final NLO result (51) 



where U{z\i) = (U(zi,z 2 )). We have calculated such 
integrals by using the representation of the type 



V? = 



'12 



2 2 



(2k • Co) 2 



based on the decomposition ( 35 1 of transverse 5- 



is defined as a coefficient function in front of composite functions. (For this example r(UJ = B( 7 , 7)r(H 



operator (17) defined with a counterterm which restores 
the conformal invariance in M = 4 amplitudes and in 
our case leads to the conformal impact factor (since the 
impact factor is given by tree diagrams it should be con- 
formally invariant even in QCD). As a consequence, the 
impact factor depends on a new parameter a (an ana- 
log of the factorization scale [i in usual OPE) which we 
chose in such a way that all the energy dependence is 
shifted in to the matrix element, leaving the impact factor 
energy-scale invariant. To compare with the result of Ref. 
[2"4"j representing the coefficient function of a usual dipole 
(without counterterm subtraction) we should trace one 
step back and look at the impact factor Iffi (z±, z 2 ,z^\ rf) 



given by ( 16 ). One should then perform Fourier transfor- 



mation to momentum space with respect to the positions 
x and y of the two electromagnetic currents in formula 
([7]) and compare it to the result (58) from Ref. [24] in- 
tegrated over z\ (and z 2 when appropriate). Hopefully, 
after these integrations the two results will coincide. 
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06OR23177 under which the Jefferson Science Associates, 
LLC operate the Thomas Jefferson National Accelerator 
Facility, and by the grant DE-AC02-05CH11231. 



VIII. APPENDIX A 

There is a subtle point in the Fourier transformation 
of the Eq. ([38]) in the forward case (cf. Ref. [20] )■ To 
illustrate it, consider the simplest term in the r.h.s. of 
Eq. ((23} 



dn a dn? 2 



dx^ dy 
2K 2 



(x - y) 



2 {dliv " 



K 2 



2 (x-yUx-y)^ (7g) 



(x - yf 



The corresponding contribution to the T-product of cur- 
rents (21) is proportional to 



/ 



dzxdzi , 
— 2 — U{z\,zi) 

Z 12 



n 2 



(x - yf 



(x-y)^{x-y) v 



(x - yf 



(79) 



7)r(2 — 7)). The integral over z\ and z 2 in Eq. (80) 
is of the form 



dzldz \Jkrf +iV M.) (81) 
-ft- (*2»)*-*7(m) 



(/// f (/: ;(/; j / z{~ • '-+"' 

2-K.j zfo Vz? n z 



where /(/i) = Jdzf 2 (z 2 2 )~i +lfl f(zf 2 ). To calculate the 
integral in the r.h.s. this equation, we take the orhtogo- 
nality condition for conformal eigenf unctions 20J 



d 2 z\ d 2 z 2 



-12 



/ z\ 2 \ g-Hiy z \ 2 x 



7Tij(5 — l/I) 



5(v - n) 

2J 1 



(82) 



and perform the inversion Zi — > %. We obtain 
'd 2 Zl d 2 z 2 / z?„ ^ r 



-12 



( Z \2 \ 5 1 2 \i-iu n 

2i[iB(\ + in) 



'10^20 ' 
1 



8(v - n) 
(83) 



so Eq. (81) turns into 

' ' dz x dz 2 / z\ 2 



-12 



{■fiirf 1 A) 

vz 10 z 20 y 



(84) 



7T 

4^2 



z ' 2vB(± + iv) 



Substituting this equa tion with f(z 2 2 ) = U(zi 2 ) 
z 2 2 V(zi 2 ) into Eq. (80) we get 



7r / a a e n-(x-v) 



. jd xdy5{y.) {x _ y)6 



9n» 



(x-y) IM (x-y) l 



(x - yf 



(x - y) 2 x*y* 



^ ^ L (x* — y* f 
B (l) \( x - y) 2 x*y* 



B(j) L (at, - y*) 5 



V{-v)) 



(85) 
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where V{v) = ^d 2 z{z 2 )-^ +w V{z). Here the first 
term comes from z$ — oo while the second from finite 
z . It is easy to see that the two terms coincide after 
change of integration variable v o — v so effectively the 
contribution of the integral over finite zq is doubled: 



J d 4 xd 4 y 5{y.) e tq < x - y ^ J 
(x - y)^ - y) v 



(x-yf 



- I d 4 xd 4 y 8(y.) 



9^ 



(x - y) e 

U(z 12 ) 
n (x-y)^{x-y) v 



{x-yf 



x r(v) 



iq-(x-y) 



B{l) {x - y) 6 L (x* - t/*) 



(x - y) 2 x*y* 



V{-v) (86) 



Since we have not used the explicit form of the Lorentz 
structure in fi and v indices, it is clear that the doubling 
effect is general for any contribution to forward Fourier 
transform of Eq. (46) (see also the discussion of zero 
transfer momentum limit in Ref. |20) ) . 
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